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Abstract 

We have calculated inclusive two-jet cross sections in next-to-leading order QCD 
for low Q 2 ep collisions superimposing direct and resolved contributions. Infrared and 
collinear singularities in the virtual and real contributions are cancelled with the phase 
space slicing method. Various inclusive two-jet distributions have been computed. The 
results are compared with recent data from the ZEUS collaboration at HERA. 



*Now at Deutsches Elektronen-Synchrotron (DESY), NotkestraBe 85, D-22607 Hamburg, Germany. 

t Supported by Bundesministerium fur Forschung und Technologie, Bonn, Germany under Contract 
05 7HH92P(0) and EEC Program "Human Capital and Mobility" through Network "Physics at High Energy 
Colliders" under Contract CHRX-CT93-0357 (DG12 COMA) 



1 Introduction 



The analysis of large transverse energy (Et) jets produced in ep scattering with very small 
virtuality (Q 2 ~ 0) offers the possibility to learn about the partonic structure of real photons. 
The hard momentum scale Et allows the application of perturbative QCD to predict cross 
sections for the production of two or more high-.&r jets. 

As is well known, in lowest order (LO) QCD two distinct processes contribute to the hard 
scattering cross sections. In the direct process, the incoming photon interacts in a point-like 
fashion with the quark from the proton via the QCD Compton scattering ^q —>■ gq or via the 
photon-gluon fusion process 7$ — > qq, if a gluon comes out of the proton. In the resolved 
process, the photon acts as a source of partons, which can scatter off the partons in the proton. 
Therefore, the cross section for the resolved process depends on the parton distributions in the 
photon in addition to the parton distribution functions in the proton. 

The two processes have distinctly different event structures. Whereas the direct process 
results in events with a characteristic (2+l)-jet structure, i.e. two high-Er jets and one Iow-Et 
remnant jet from the proton, the resolved photon events have a characteristic (2+2)-jet struc- 
ture, where in addition to the two high-.&r jets and the low- Et proton remnant a second jet of 
low-.&r fragments of the photon is produced. While in the (2+l)-jet events the total photon 
energy contributes to the hard scattering process, only a fraction participates in the (2+2)-jet 
events. 

This characterization of the direct and resolved processes is valid only in LO. In next- 
to- leading order (NLO), the direct cross section may also have contributions with a photon 
remnant jet. In addition, both components are related to each other through the factoriza- 
tion scheme and scale at the photon leg. The dependence of the NLO direct cross section on 
this scheme/scale must cancel up to NNLO terms against the scheme/scale dependence of the 
resolved cross section. Therefore in NLO both components must be considered together and 
consistently be calculated up to NLO in the parton distributions of the photon and the proton, 
respectively, and in the hard scattering cross sections. 

The simplest observable is the inclusive single-jet cross section. It depends on the transverse 
energy and the rapidity of the observed jet. Complete NLO calculations for this cross section 
have been done previously J]], |2], |3|] and have been compared to experimental data from ZEUS 
|U and HI [|]. In these calculations special techniques for cancelling soft and collinear singular- 
ities have been applied, which are suitable only for the computation of inclusive single-jet cross 
sections. To obtain different cross sections, as for example the inclusive two-jet cross section, 
the whole calculation must be repeated. Inclusive two-jet cross sections depend on one more 
variable. Therefore they are a much more stringent test of the QCD predictions than inclusive 
single-jet cross sections. Furthermore, for a comparison with experimental data the calculations 
must be performed in such a way that the experimental cuts can be built in easily. For this 
purpose, the so-called phase space slicing method is suited particularly well. In this method, 
an invariant mass resolution cut is introduced to isolate the soft and collinear singularities 
of initial and final states. This resolution cut is purely technical and is chosen small enough 
not to disturb any experimental constraints or jet recombination requirements. Based on this 
method we recently calculated the inclusive two-jet cross section for direct photoproduction ||. 
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The results of this calculation were compared with data from the ZEUS collaboration J7|, in 
which the direct contribution and the resolved contribution were separated with a cut on x 7 , 
i.e. the momentum fraction of the photon entering the hard scattering. The contribution of 
the resolved photon in the enriched direct 7 sample was estimated in lowest order since it was 
supposed to contribute only a fraction in this sample. As a result of this investigation it turned 
out that the resolved contribution in the enriched direct 7 sample was not negligible and could 
be, depending on the assumed photon structure function, as large as 30% of the total two-jet 
cross section. It is clear that the estimate of the resolved cross section in this region by a LO 
computation is insufficient and that the NLO corrections of the resolved cross section must be 
included. In addition, in the same analysis of the ZEUS two-jet data |7[ the inclusive two-jet 
cross section in an enriched resolved 7 sample was measured, which could be a decisive test 
of the available parton distributions of the photon. For this reason and since we expect more 
two-jet experimental data in the future, a complete NLO calculation of the inclusive two-jet 
cross section is mandatory. For the direct component, we rely on our previous work |6], ||. 
The NLO calculation of the resolved component has been completed recently |§. A detailed 
description of the calculation will be presented in a separate publication JED]. Here we shall 
only give a short presentation, report some results on various inclusive two-jet cross sections, 
which might be measured in the future, and recalculate the two-jet cross section in the enriched 



direct 7 sample and the enriched resolved 7 sample for a comparison with recent ZEUS data [11 



The outline of the paper is as follows: in section 2, we explain the formalism used to 
calculate the inclusive two-jet cross section for the resolved part. In section 3, we present some 
results in order to demonstrate the flexibility of our method and a comparison with recent 
ZEUS measurements. Section 4 contains a summary. 



2 Next-To-Leading Order Cross Sections 
2.1 Photon Spectrum 

Before we explain how the NLO corrections have been calculated, we must specify the relation 
between the ep and the jp cross sections. Let us start with the electroproduction process 

e(k)+p(p)^e'(k')+X, (1) 

where k, k', and p are the four-momenta of the incoming and outgoing electron (positron) and 
the proton, respectively. X denotes a generic hadronic system which will be specified later. 
q = k — k' is the momentum transfer of the electron to the photon with virtuality Q 2 = —q 2 ~ 0. 
For small Q 2 , the ep cross section can be calculated in the Weizsacker- Williams approximation, 
where the 'jp cross section can be factorized and the photon spectrum is given by the function 



a 1 + (1 - x a ) 2 Q 2 mSLX (l - x a) 

fi-y/eKXa) ~ ~ Jn — . (2J 

Z7T X a m e X a 



Here, m e is the electron mass, and x a = (pq) / (pk) ~ E 1 / E e is the fraction of the initial electron 
energy transferred to the photon. The cross section for the process (1) ep — > e'X is then given 
by the convolution 

& a , max 

da(ep e'X) = J dx a F 7/e (x a )dcr(7p -> X). (3) 
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Qmax) x a,max < 1, and x ajmin are determined by the tagging conditions of the experiments at 
HERA and will be specified when we present our numerical results. In (3), da^p — > X) is 
the photoproduction cross section with initial photon energy E 1 = x a E e and photon virtuality 
Q 2 ~ 0. 



2.2 Jet Cross Sections in NLO 

As we explained in the introduction, inclusive jet cross sections da^p — > jets) receive contribu- 
tions from two components, direct and resolved, depending on the way in which the incoming 
photons participate in the hard scattering subprocesses. In the resolved case, the NLO inclusive 
two-jet cross section is calculated from the formula 



3 

]T j dy a J dx b F ah {y a ,M 2 )F b/p {x b ,M 2 p ) 

K ab {R,M^M p ^)\. (4) 



dV 



d#rd77id7/2 a , b=q , g 

( fd 3 a(ab — > jet x + jet 2 ) \ ol s (h) 



dE T dr]id7] 2 /t 2tt 



In (4), Et is the transverse energy of the measured or trigger jet with rapidity rjx. r]2 denotes 
the rapidity of another jet such that in the three-jet sample these two jets have the largest 
and second largest Et, i.e. Et 1 ,Et 2 > Et 3 - The first term on the right-hand side of (4) is 
the LO cross section, and the second term is the NLO correction. The latter depends on the 
factorization scales M 7 and M p at the photon and proton vertex, at which the initial state 
singularities are absorbed into the parton distribution functions of the photon and the proton. 
fj, is the renormalization scale. The variable R is the usual cone size parameter, which defines 
the size of the jets with transverse energy Et, rapidity T], and azimuthal angle 0. When two 



partons fulfill the Snowmass condition |12J with the cone size parameter R, they are combined 
into one jet. The same jet definition must be used in the analysis of the experimental data, if a 
meaningful comparison between theory and experiment is intended. It is clear that either of the 
two jets can consist of two partons inside a cone with radius R. Essentially the same formula 
applies for the direct cross section. In this case, the photon structure function F a /^(y a ,M 2 ) is 
replaced by 5(1 — y a ). 



The NLO corrections K ab in (4) for the hard scattering subprocesses are calculated with the 
phase space slicing method using an invariant mass cut for the separation of the cross section 
for 7p — > 2 jets and ^p — > 3 jets. This invariant mass cut is defined by IpiPj < ys, where 
s is the partonic center-of-mass energy squared. So, for example, in the resolved subprocess 
QiQj ~^ QiQjQ^ where qj can be identical or non-identical quarks, the cross section has soft, 
initial, and final state collinear singularities. The cross section da(qiq~j — > qiq^g) is integrated 
over these singular regions up to the cut-off y, which isolates the respective singularities in 
terms of poles in e = (4 — d)/2, where d is the dimension in the dimensional regularization 
method. These singularities cancel against the singular contributions of the virtual corrections 
to qiq~j — > qtljj and against the subtraction terms at the scales M 7 and M p , which are absorbed 
into the parton distribution functions of the photon and the proton, respectively. Outside the 
cut-off region, i.e. for 2piPj > ys, we have genuine qiq^g final states. The cross section in this 
region is subdivided into the two-jet cross section and the three-jet cross section contribution, 
depending on whether two of the final state partons are combined according to the Snowmass 
condition or not. The LO contribution, the NLO virtual contributions, and the NLO correc- 
tions inside the y-cvX contribute to the two-jet cross section together with the contributions 
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inside the recombination cone with radius R. The part of the qiq^g cross section not fulfilling 
the cone recombination condition is the three-jet cross section, from which we have calculated 
the inclusive two-jet cross section as a function of E T , rjx, and r] 2 . For the exclusive two-jet 
cross section, E? x = Et 2 = Et- To separate double singularities in the qiq^g cross section part 
we have used the method of partial fractioning. The same partial fractioned expressions are 
applied for the calculation inside and outside the y-cnt. Inside the y-cut region, the integrations 
are done analytically with the approximation that terms 0(y) are neglected. This is necessary 
since otherwise with a d- dimensional space-time the necessary integration could not be carried 
out. Outside the y-cut region, the contributions to the inclusive cross section are evaluated 
with no further approximations. Because of the approximation in the analytic part inside the 
y-cut, the parameter y must be chosen very small. It turned out that for y of the order 10 -3 
to 10~ 4 the inclusive cross section is independent of the cut-off parameter y. However, the 
inclusive two-jet cross section depends on the cone size R, which must be chosen in accordance 
with the analysis of the experimental data. In the case of all the other hard subprocesses for the 
resolved cross section, of which there are many and which will not be listed here, the calculation 
proceeds in the same way. The complete list will be given in a separate publication [It]] and 



can also be found in 0, [Tj|. In these references, the details of the analytic integrations and the 
cancellation of the singularities for NLO corrections is described as well. For the direct cross 
section, the calculation of the NLO correction was already reported and is given with all details 
in our previous work [|, |[ [J. 

Before the final results were obtained, some tests of the NLO corrections have been per- 
formed. For similar tests of the NLO corrections of the direct cross section see 0. The same 
tests were done for the resolved subprocesses. So we checked that the cross sections are in- 
dependent of the cut-off y if y is chosen small enough. This was the case for y < in 
all considered cases under different kinematical conditions. For y > 10~ 3 , we observed some 
small y dependence caused by the approximations in the analytical integrations. Second, the 
resolved inclusive one-jet cross section was calculated and compared with earlier results for 
which a completely different method for cancelling the infrared and collinear singularities was 
applied 0, Very good agreement was found in all channels and for the total resolved cross 
section. Third we tested that the sum of the NLO direct and the LO resolved cross section is 
independent of the factorization scale M 7 similar to such a test for the single-jet inclusive cross 



section earlier 15 . Details of these tests can be found in 



For the calculation of the inclusive jet cross sections, we need the parton distributions F a / 7 
of the photon and Fb/ P of the proton, respectively (see (4)). We have chosen the NLO set of 
Gliick, Reya, and Vogt (GRV) in the MS scheme for F a / 7 [[Uj and the CTEQ3M set, i.e. also 
a MS set in NLO, for Fb/ P Hl7| . All scales are set to /i = M 7 = M p = Et- ot s is calculated 
from the two-loop formula with rif = 4 massless flavors with = 239 MeV, equal to the 
A-value of the CTEQ3M parton distributions. The charm quarks are treated as light flavors 
with the boundary condition that the charm content of the two structure functions vanishes 
for Ml, Ml < m 2 c with m c = 1.5 GeV. 
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3 Results and Comparison with ZEUS Data 



In this section we present some characteristic numerical results for two-jet cross sections con- 
centrating on the resolved contribution. All results are obtained with the phase space slicing 
method. Except where we compare with recent ZEUS data, we do not apply special cuts on 
kinematical variables of the initial or final states dictated by the experimental analysis, although 
this is necessary if one wants to confront the prediction with experimental data. Therefore, 
the photon energy fraction x a in (2) is allowed to vary in the kinematically allowed range from 
x a ,min < x a < 1 • <5max = 4 GeV 2 in (2) is the same as usually used in the jet analysis of the 
ZEUS experiment. 

The further evaluation of the two-jet cross section is based on two separate contributions - 
a set of two-body contributions and a set of three body contributions. Each set is completely 
finite, as all singularities have been cancelled or absorbed into structure functions in the same 
way as in the case of the direct cross section ||. Each part depends separately on the cut-off 
y. The analytic contributions are valid only for very small y, since terms of 0(y) have been 
neglected. For very small y, the two contributions have no physical meaning. In this case, 
the terms depending on In y force the two-body contributions to become negative, whereas the 
three-body cross sections are large and positive. The resolved cross section shows a similar be- 
havior || . When the two- and three-body contributions are superimposed to yield the inclusive 
two-jet cross section, the dependence on the cut-off will cancel. This will only be the case if 
the inclusive cross sections are defined in the correct way, so that they are infrared safe. We 
have checked this explicitly by varying y in the superposition of the two pieces for the inclusive 
one-jet and two-jet cross section separately. For the two-jet cross section, it is essential that the 
variable Et is chosen in the way as stated above, namely that Et is the transverse momentum 
of the two largest Et jets. This means that the cut-off y is purely technical. It only serves to 
distinguish the phase space regions, where the integrations are done analytically with arbitrary 
dimensions from those where they are done numerically in four dimensions. When one com- 
pares to experimental data with special kinematical cuts, y must be chosen sufficiently small 
in order to make sure that the cuts on kinematical variables of the final state do not interfere 
with the cancellation of the y dependence. 

For the recombination of two partons i and j with the cone constraint Rij < R, where 



and rjj, <pj are the rapidity and azimuthal angle of the recombined jet, we choose the radius 
R — 1. In some cases, an ambiguity may occur when two partons i and j qualify both as two 
individual jets i and j and as a combined jet. In this case we only count the combined jet to 
avoid possible double counting. It is clear that in NLO the final state may consist of two or 
three jets. The three-jet sample consists of all three-body contributions, which do not fulfill 
the cone constraint given above. 

First we show results for the resolved inclusive two-jet cross section in fig. [I]. We fixed rji 
at r)i = 1 and 772 at three different values of 772 = 0, 1, 2 and plot the cross section as a function 
of Et = Et x - The curves for 7/ 2 = and 772 = 2 are rescaled by factors of 0.1 and 0.5. Since 
the two-jet cross section is much more exclusive than the one-jet cross section, it is reduced in 
magnitude as we expect it. In fig. 1, we also show the LO cross section as a function of Er y . 




(5) 
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Figure 1: Inclusive dijet cross section d z a j dET 1 dr\\dr\ 2 for resolved photons as a function of 
for rji = 1 and three values of r\ 2 = 0, 1,2. The cross section for i] 2 = (r] 2 = 2) is multiplied 
by 0.1 (0.5). 




Figure 2: Inclusive dijet cross section d 3 a/ dE Tl dr]idr] 2 at E Tl = 20 GeV for resolved photons, 
as a function of r/i and r\ 2 . The LO plot (left) is exactly symmetric, the NLO plot (right) only 
approximately. 
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The ratio of NLO to LO is 1.7 over the whole E? l range. In this comparison, the LO cross 
section is calculated with the LO subprocesses, but with the same definition of a s and NLO 
structure functions. For r/ 2 = the NLO and LO differ very much for large E^- In this case 
the third jet in NLO opens up additional phase space, so that Et x can go up to Et x < 50 GeV. 
A similar behavior is observed for the corresponding direct cross section shown in ||. 

Next we present the dependence of the resolved cross section on the two rapidities in form 
of the three-dimensional lego-plots in fig. |2|. The equivalent plot for the direct cross section was 
already given in |]]. The leading order is shown on the left side and is completely symmetric 
in 771 and r/ 2 . The NLO cross section, which is shown on the right-hand side of fig. |^, is not 
symmetric any more due to the presence of a trigger jet with transverse energy E^, which is 
fixed at Et-l = 20 GeV. This can best be seen at the bottom of the contour plots, where at least 
one of the two observed jets is far off the central region. The NLO predictions are considerably 
larger than the LO predictions. For the direct cross section, the NLO cross section is only 
moderately larger than the LO cross section |§, where the fc-factor is 1.25. 

This becomes even clearer when we plot the projections of the lego-plots for fixed r/i — 0, 1, 2 
and 3. In fig. |^, the LO and NLO distributions in r/ 2 are shown for the direct photoproduction 
and in fig. ^ for the resolved part. In the direct case it is clearly seen that the second jet tends 
to be back-to-back with the first jet, since the maximum always occurs at r/2 — r/i. However, 
at 771 = 3 this is no longer possible due to phase space restrictions. The 772-distributions for 
resolved photons in fig. |4] are considerably broader than in fig. ||] due to the smearing of the 
hard cross sections with the distribution functions of the partons in the photon. The maxima 
of the curves are also not so much dominated by kinematics but more by the quark and gluon 
structure of the photon in different x 1 regimes. They do not lie at r/2 = f]i- Therefore, two-jet 
rapidity distributions are very well suited to constrain the photon structure. The fc-factors 
range from 1.65 to more than 3 in the very forward region of the proton. The shapes of the 
distributions are very similar in LO and NLO. The absolute values make, however, an impor- 
tant difference. 

Next we show the sum of the direct and resolved cross sections which gives the physical 
complete photoproduction result. First we have plotted the two-jet cross section as a function 
of the transverse energy E Tl at 771 = 772 = 1- Fig. || gives the LO result, fig. |6| the NLO result. 
As we can see, the point where direct and resolved contributions are equally important is near 
Et x = 20 GeV in leading order and Et x = 30 GeV in next-to-leading order. These crossing 
points are somewhat smaller than for the inclusive one-jet cross section at 77 = 1, so that direct 
photons make a stronger impact in dijet production. 

If one plots the complete two-jet cross sections as a function of r/ 2 , the different behav- 
iors of direct and resolved photons add up to the full curves in fig. ^ (LO) and fig. || (NLO). 
These plots are best suited to decide in which rapidity regions one can see best the resolved 
photon structure. We have already seen that this will be in situations where the two jets are 
not back-to-back, e.g. for r/i = and positive r/ 2 in the upper left plots of fig. ^ and |8|. On 
the other hand, the proton structure can best be studied with direct photons, when the cross 
section is not folded with another distribution. A possibility is at 771 = and negative values 
of 772. This is especially interesting for the small- a; components of the proton like the gluon 
and the quark sea. Another interesting observation is that the relative importance of direct 
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Figure 3: Projections of the LO (full curves) and NLO (dashed curves) triple differential dijet 
cross section for direct photons at E Tl = 20 GeV and fixed values 0/771 = 0, 1, 2, and 3, as a 
function of r) 2 . 




Figure 4: Projections of the LO (full curves) and NLO (dashed curves) triple differential dijet 
cross section for resolved photons at E Tl = 20 GeV and fixed values of i]i = 0, 1, 2, and 3, as 
a function of r] 2 . 
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Figure 5: Inclusive dijet cross section d 3 a/ dE Tl dr]idr] 2 for full photoproduction at r\\ — rj 2 — 1 
as a function of E Tl . The full curve is the sum of the LO direct (dotted) and LO resolved 
(dashed) contributions. 
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Figure 6: Inclusive dijet cross section d 3 a/ dE Tl dr)idr)2 for full photoproduction at rji = r] 2 = 1 
as a function of E Tl . The full curve is the sum of the NLO direct (dotted) and NLO resolved 
(dashed) contributions. 
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Figure 7: Projections of the complete triple differential dijet cross section at E Tl = 20 GeV and 
fixed values 0/771 = 0, 1, 2, and 3, as a function of r/ 2 . The full curve is the sum of the LO 
direct (dotted) and LO resolved (dashed) contributions. 




Figure 8: Projections of the complete triple differential dijet cross section at E Tl = 20 GeV and 
fixed values 0/771 = 0, 1, 2, and 3, as a function of r} 2 . The full curve is the sum of the NLO 
direct (dotted) and NLO resolved (dashed) contributions. 
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and resolved processes changes when calculating dijet photoproduction in next-to-leading order 
0{aa 2 s ): resolved processes are much more important at E Tl = 20 GeV than one would have 
guessed from a leading order estimate. 



For the comparison of our two-jet predictions with experimental data it is essential that for 
the theoretical calculations the same jet definitions are introduced as used in the experimental 
analysis. First experimental data for inclusive two-jet production have been published by 
the ZEUS collaboration in and [pl|. The analysis of these data was based on the usual 



cone algorithm with R = 1 [^] similar to the cone algorithm which we used in our NLO 
calculations presented above. The recent ZEUS analysis based on the 1994 data taking 
extends the earlier analysis based on the 1993 data in several ways. The larger integrated 
luminosity utilized in 1994 lead to a reduction of the statistical errors as well as allowing 
for the measurement of the cross section at higher Et, a region, where uncertainties due to 
non-perturbative hadronization of partons into jets are expected to be reduced making the 
comparison with the NLO predictions more meaningful. Furthermore, the ZEUS collaboration 
applied three different jet definitions: two variations of the cone algorithm |T2[ are used which 
they called "EUCELL" and "PUCELL" . These algorithms treat seed finding and jet merging 
in different ways. In addition, the kj< cluster algorithm "KTCLUS" W§ 
collisions is used. In 



for hadron-hadron 



IE 



several versions of the kr algorithm have been introduced which 
left some flexibility about how the stopping condition and the recombination scheme can be 
implemented. A version which is particularly suitable to define jets for inclusive jet cross section 
measurements was emphasized in f20| . It is also closest to the cone algorithms and states that 
two protojets i and j (i.e. two partons in the three parton final states) are merged if 



\l (7H ~ Vj] 



< 1 



(6) 



and the recombination of the two protojets is done with the E? scheme as is usually also done 
in the cone algorithms [12 1. Compared to the cone algorithm [1^] with R — 1, which was 
applied so far to generate the results in figs. 1-8, the two partons are merged if they are less 
than R = 1 apart in rj — space, whereas in the cone algorithm the merging condition is on 
the opening angle between either of the two partons and the jet center obtained from the Et 
recombination equations, i.e. Ri < R and Rj < R with R = 1 and Ri defined in equation 
(5). Thus this version of the kx algorithm is very similar to the cone algorithm usually applied 
in hadron-hadron and photon-hadron collision. As it is a cluster algorithm, the ambiguities 
associated with seed finding and jet merging are avoided. The same k T algorithm, denoted by 



KTCLUS, was used by the ZEUS collaboration in the analysis of their recent data [IT|. 



Except for the different jet definition, the cross section to be compared to the ZEUS data is 
the same as in our earlier work ||. It is the inclusive two-jet cross section d 3 a/dETdr]dr]*, where 
rj — |(?7i + 772) is the average rapidity of the two observed jets and Et is the transverse energy 
of the "trigger" jet. This cross section is integrated over the difference in rapidity, \rf\ < 0.5, 
and integrated over Et 1) Et 2 > E T im with varying E T im . The additional restrictions on the 
individual rapidities rji and r^ 2 used in the ZEUS analysis have not been implemented since 
they only modify the cross sections at r\ ~ 2. As the experimental constraint on the transverse 
energies of both jets is not infrared safe in NLO, we allow the second jet to have a transverse 
energy less than E^™ if the third unobserved jet is soft, i.e. has a transverse energy of less than 
1 GeV (Et 3 < 1 GeV). Through this procedure we avoid the dependence of the theoretical 
prediction on the y-cut. In addition, we separate "direct" and "resolved" contributions with 



11 





Figure 9: Inclusive dijet cross section dcr/drj at x® BS > 0.75 as a function ofrj and integrated 
overt]* G [—0.5,0.5] and Ex > 6, 8, 11, and 15 GeV. Our leading and next-to-leading order 
predictions are compared to preliminary 1994 data from ZEUS using the KTCLUS algorithm. 





Figure 10: Inclusive dijet cross section da/drj at x® BS G [0.3,0.75] as a function of Tj and 
integrated overt]* G [—0.5,0.5] and Et > 6, 8, 11, and 15 GeV. Our leading and next-to- 
leading order predictions are compared to preliminary 1994 data from ZEUS using the KTCLUS 
algorithm. 
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the variable 

obs Ei^e-w 



7 ~ 2x a E e ' ^ 



which measures the fraction of the photon energy that goes into the production of the two 
hardest jets. As in the ZEUS analysis the enriched direct 7 sample is defined with the cut on 
x obs > 0.75. The enriched resolved 7 sample is obtained for 0.3 < x^ BS < 0.75. The very 
low x® BS are excluded in the experimental analysis since this region is not accounted for by 
the Monte Carlo routines necessary to correct for detector effects. Fig. [| shows the "direct" 
cross sections d<r/dr/ as a function of fj in LO and in NLO for E^ in = 6, 8, 11, and 15 GeV 



together with the preliminary ZEUS data (TT|. The LO and NLO predictions are very similar. 
We emphasize that the LO curve is calculated with the same NLO structure functions as in 
the NLO calculations, only the hard scattering parton-parton cross sections are evaluated in 
LO. The NLO curves are shifted to larger rj and agree in average better with the data than 
the LO curves. The general agreement is quite good indicating that the GRV choice for the 
photon parton distributions needed for the resolved contribution is quite reasonable in the high 
x 7 region. In our previous publication 0, we compared the corresponding 1993 ZEUS data to 
our prediction, which then only included the LO resolved contributions for the GRV photon 
parton densities. 

The "resolved" photon two-jet cross section, where 0.3 < x^ BS < 0.75, is compared to the 
ZEUS data in fig. [ITj for the same E™ m values. Here the agreement is satisfactory only with the 
data in the high E? regions E™ m = 11 and 15 GeV. At the lower Et, the experimental cross 
section is larger than the NLO prediction. The disagreement increases with decreasing E^ n . 
We attribute this difference between theory and experimental data to additional contributions 
due to multiple interactions with the remnant jet not accounted for by our NLO predictions. 
As one would expect that these effects diminish with increasing Et, the agreement between 
theory and experiment improves in the large Et region. We also studied the uncertainty com- 
ing from the insufficiently constrained gluon in the photon by recalculating the resolved cross 
section with twice the GRV gluon distribution. At low Et, the cross section is enhanced by 
33% thus reducing the discrepancy, whereas at large Et it is only enhanced by 20%. With 
more accurate data at Et > 17 GeV, such a gluon distribution could be excluded. In addi- 
tion to the experimental errors shown, which include systematic and statistical errors added 
in quadrature, there is a systematic uncertainty arising from the uncertainty in the calorime- 
ter energy scale, which is highly correlated between bins and is therefore excluded from the 
systematic errors shown in figs. and [TO. This uncertainty is largest for the "resolved" cross 



section at i?™ m = 6 GeV and leads to a ±0.5 nb uncertainty in the cross section near 77 ~ 0.6 



rjj. In the "direct" cross section, this energy scale uncertainty is smaller by a factor of two |TTJ . 



In fig. pi we observe that the NLO two-jet cross sections are somewhat smaller than the LO 



cross sections, which are included only to see the effect of the NLO corrections. These LO cross 
sections are independent of the jet definition, i.e. are always the same for the cluster algorithm 
or a cone algorithm with any cone radius. They should not be compared to the experimental 
data since it is known that the two-jet cross sections depend on the jet definitions. The choice 
of the jet definition has an effect of about 25-30% in both data and theory in the "direct" cross 
section and of about 50% for the "resolved" cross section [lj, [2l| . 
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4 Summary 



Differential cross sections d 3 0"/d-Erd?7id?72 have been calculated in NLO for direct and resolved 
photoproduction. Infrared and collinear singularities are cancelled with the phase space slicing 
method using an invariant mass cut-off. With this method we are able to incorporate various 
cuts on the final state as used in the analysis of experimental data and to perform calculations 
for different choices of jet algorithms. Numerical results for the two-jet inclusive cross sections 
at HERA have been presented employing the usual Snowmass cone algorithm for the jet def- 
inition. For a cone radius of R = 1, the NLO corrections lead to an increase of the order of 
70% compared to the LO prediction in the resolved contribution and of the order of 20% in the 
direct case already presented earlier M. 



Using a particular version of the kx cluster algorithm [2~0j, we calculated the cross section 
da/drj with cuts on x° BS to separate "direct" and "resolved" contributions as in the analysis 



of the ZEUS data J7|, |TT|, |TB|. We find fairly good agreement with recent ZEUS data [[□] 
for the enriched direct 7 sample. For the enriched resolved 7 sample, the agreement is good 
in the larger Et region, _E™ m > 11 GeV. For smaller Ej- (E™ m = 6 and 8 GeV), we find a 



discrepancy which we attribute to additional multi-parton interactions [|ll] not incorporated in 
the theoretical calculations. 
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